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Abstract 

Theoretical cross sections for the pressure broadening by hydrogen of rotational transitions of 
water are compared to the latest available measurements in the temperature range 65—220 K. A 
high accuracy interaction potential is employed in a full close coupling calculation. A good agree- 
ment with experiment is observed above ~80 K while the sharp drop observed experimentally at 
lower temperatures is not predicted by our calculations. Possible explanations for this discrepancy 
include the failure of the impact approximation and the possible role of ortho-to-para conversion 
of H 2 . 



Water is a prominent molecular component of interstellar matter. It has been observed in 
most astrophysical environments, both in gas and solid phase, as the most abundant poly- 
atomic molecule [TJ. Understanding water spectra is a key to the general thermodynamical 
budget of astrophysical objects, because of the many allowed spectral transitions, in mm, 
sub-mm and infra-red regions of the spectrum. Also, the chemical and even exo-biological 
importance of water cannot be understated. A detailed comprehension of water physical 
chemistry in various environments is a major goal of the newly launched Herschel Space 
Observatory (HSO). 

In order to extract information from a molecular rotational spectrum, it is essential to 
model its excitation scheme. Indeed, at low densities, some lines may appear in absorption, 
some others in emission and there is no a priori reason for the molecule under scrutiny to be 
in thermodynamical equilibrium with the main neutral gas, H2. Hence retrieving physical 
information from spectral lines is only possible with a careful modeling of the interaction of 
the water molecule with its photonic and molecular environments. Obtaining such models 
has been a continuous effort in three directions: studies of radiative transfer mechanisms, 
scattering theory, and very importantly, laboratory experiments capable of testing various 
theories and models pertaining to collisions [2]. Many comparisons between experiments 
and theory are nowadays underway, in order to put the water-hydrogen interaction on a 
firm ground [3H6] . 

In this communication, we wish to show the first results of a fully ab initio computation 
of pressure broadening cross sections, o" p b '(T), based on a high precision potential energy 
surface (PES) for the water-hydrogen interaction [3 |8] (hereafter V08), as compared to the 
very recent experiments of pressure broadening at low temperatures by Dick et al. [9j 
[TO] (hereafter DDP10). To our knowledge, it is the first time that such a full ab initio pressure 
broadening quantum calculation is performed for a non-linear molecule in interaction with 
another molecule. Other comparisons were very successful for simpler symmetries [TT1TI5] . 

Water being an asymmetric rotor, the rotational levels are usually denoted as j KaKc or 
(j, r) where j is the rotational quantum number associated with the angular momentum and 
(/€ a , k c ) (projections of j along inertia axis) are pseudoquantum numbers, and r = K a — k c . 
The rotational constants of H 2 are taken at A = 27.88063134 cm" 1 , B = 14.52176959 cm" 1 , 
C = 9.277708381 cm -1 . Rotational constant of H 2 is taken at B = 60.853 cm" 1 . 

Following the experiments of DDP10, we compute here the pressure broadening of the 
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two spectral transitions connecting the ground states of water: the para 1113 GHz line 
(In o) and the ortho 556 GHz line (li ^— loi)- All our calculations are based on 
the V08 water-hydrogen full-dimensional PES which was obtained by combining standard 
CCSD(T) calculations with elaborate explicitly correlated CCSD(T)-R12 calculations. As 
in Dubernet et al. [16J, we have employed the rigid-body version of the V08 PES obtained 
by averaging the full-dimensional PES over the ground vibrational states of the monomers. 
Full details can be found in [8]. 

The broadening of a rotational spectral line because of collision with a buffer gas has 
been studied theoretically and experimentally for a long time and theory is by now well 
established. The very general impact approximation states that collision times are much 
shorter than time between collisions. Within that approximation, which we discuss later, 
it has been shown that the pressure broadening cross section for the transition from initial 
state i to final state /, at temperature T, cr^'j(T), may be expressed by closed expressions 
based on the transition matrix T [17, 18J. 

For a broadening coefficient Tj^iiT), in frequency /pressure units, Baranger [Tj5] defined 
the pressure broadening cross section at energy E, a p ' h \E) as 

IWT) = 1/2 (nvaf^E)) (1) 

1 V ' 1 \ IT Boltzmann av. V ' 

where % and / are the initial and final states of the transition, n is the density of the 
observed molecule, and v is the relative velocity of water and hydrogen. E is the collision 
kinetic energy. Hence it is possible to define a Boltzmann averaged cr^'^T): 

<i( T ) = ^J <me~ E/T E dE. (2) 

Two equivalent ways have been proposed to calculate a^^E), and consequently the 
averaged Oy^(T'). Following Baranger [IT] . Schaefer & Monchik [201 121] and Green [22], 
we have, for a rotational transition of H 2 (Jf, Tf <— ji, Tj), assuming that H 2 remains in an 
unchanged j 2 state: 



•jr J _ _ 

<T P ' h i3fTfj2 <~ jiTij 2 ;E) = £ XfaJfJiJfjJj-uj'uJviju) 



(3) 



J i J f I,l',jl2j' 12 ,jl2j[ 2 

x (32loJi2jiri\T J ^E i )\j 2 l' f 12 j l T i )(j 2 l ] 12 j f T f T Jf (Ef) 3 2 l' f 12 3 f T f ) 
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In Eq.g), ji 2 ,j[ 2 are the angular quantum numbers numbers resulting from the coupling of 
angular momenta j 2 and j { / j f ( for example, \j 2 - j { \ < j 12 , f 12 < (j 2 + ji))- kJ'o are the 
orbital quantum numbers. T J (E) = 1 — S J (E) is the transition matrix, at total angular 
momentum J. The X(.) function groups all angular coupling coefficients and parity sign 
terms; it is explicit in [21]. Ei and Ef are the two initial and final total energies (Ei ^ Ef), 
corresponding to the same kinetic energy, (hk) 2 /2fi, /i being the collision reduced mass and k 
the momentum. A similar equation relevant for H 2 0- H 2 coupling may be found in Monchik 
[21], Eq. (1). It should be noted the l/(2j 2 + 1) factor in front of Eq.g, lacking in [22] is 
necessary, as underlined for example in |23j . 

Baranger, followed by many authors, proceeded to compute the a p ' h '(E) in a different, yet 
equivalent way, with help of the optical theorem. We have, with nearly the same notations 
as 02! : 
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E^(£) + E4 n ^(£) 

+ J \m;E)-f f (n-,E)\ 2 <m (4) 



In Eq.Q, a mel (E) are ordinary inelastic cross-sections, i' and /' being all levels connected 
to / or % at kinetic energy E. The fi(Q;E), ff(Q;E) are the differential elastic scattering 
amplitudes, for the initial and final states, which must be set to interfere before taking the 
modulus and integrating over all scattering angles Q. 

It must be strongly underlined that both approaches are equivalent. Quite often, since 
inelastic cross sections or rate coefficients are made available in the literature, in order to 
have an estimate of a^^E), equation (4) is truncated: only the first two terms, the inelastic 
cross sections, are used, yielding sometimes to reliable results [21] and sometimes not [13J, 
depending on the structure of the scattering amplitudes (see below). 

In order to calculate the T matrix elements of Eq. ([3]), we performed a full quantum 
close coupling scattering computation with help of the Molscat program [25]. The T ma- 
trix elements were subsequently summed at each kinetic energy to get the relevant o~(E) 
cross-sections, inelastic and pressure-broadening. We separately computed collisions of the 
four symmetry types: (para/ortho H 2 0)-(para/ortho H 2 ). Parameters of the Molscat cal- 
culations are similar to the ones used previously [HI [323 I2SH2B], with the following rotational 
basis sets: para-H 2 : j 2 = 0,2; ortho-H 2 : j 2 = 1. Ortho and para H 2 0: ji < 5,6,7, for 
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increasing collision energy. For E < 20 cm -1 , the hybrid modified log-derivative/ Airy prop- 
agator of Manolopoulos and Alexander was used; above that energy, the diabatic modified 
log-derivative method of Manolopoulos was used [25] . We checked convergence with respect 
to basis set size, maximum range of radial integration and size of step in the radial propaga- 
tion. The collision energy range was 0.5 < E < 350 cm -1 , with increments AE small enough 
to roughly describe resonances. Decrease of the AE step did not change significantly the 
averaged cr^'^T), which was obtained by a numerical integration of Eq. (2). At collision 
energies E > 350 cm -1 , a p ' h '(E) remains essentially flat and was therefore extrapolated as 
a constant. 

Results are presented in figures [I] and [2j In addition to the pure ortho H2 (blue lines) 
and para H 2 (red lines) cross sections, we added two possibilities for the ortho-to-para ratio 
(OPR) of H2: the solid black lines suppose a "normal" OPR value of 3, as expected in the 
DDP10 experiment. The grey solid lines suppose that the thermodynamical equilibrium 
OPR value is reached at each temperature T, as if ortho to para transitions were possible. 
In both figures, the results of DDP10 are shown as green open symbols, with their exact 
values provided courtesy of B. Drouin. 
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FIG. 1: Pressure broadening cross-sections for the transition at 556 GHz, ortho-E^O, lio loi- 
Theory: Red line, below, pure para-H2 collisions; blue line above, pure ortho-H2 collisions. Black 
continuous line, OPR value, 3; grey dashed line OPR value at equilibrium for each T. Green 
triangles, experimental values taken from [10J. 



Several points may be seen by inspection of figures [T] and |2j If we assume an OPR value of 
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FIG. 2: Pressure broadening cross sections for the transition at 1113 GHz, para-H^O, In <— Ooo- 
Same caption as for figure [TJ 

3, we see that the theory vs. experiment agreement is very good (within 30%) for T > 80 K. 
Let us recall that there are no adjustable parameters involved in the simulation, except for 
the OPR value. This agreement should come at no surprise, since recent analogous calcu- 
lations on branches of Raman spectra show similar successes [TTJ 1221 ED]- Both these 
calculations and ours make use of state-of-the-art ab initio PES as well as fully converged 
close-coupling calculations. Within the impact approximation valid at these temperatures 
and moderate densities (n ^ 1 amagat), Baranger formalism (17) . equation (J3j) , is essentially 
exact for nonoverlapping lines. The whole uncertainty that arises is due either to inade- 
quacies of the impact approximation, which is to be ruled out here at T > 80 K (see for 
example [18] or [31] for a discussion) or else to imprecision of the PES. We see thus that we 
actually test the PES by comparing experimental pressure broadening with careful quantum 
calculations, in the relevant density and temperature ranges. 

Pressure broadening cross sections, at these intermediate temperatures, are sensitive to 
the overall shape of the PES, and especially to the isotropic part. This may be understood in 
two ways. First, the T matrix elements actually used in Eq.([3]) are elastic in the rotational 
quantum numbers, thus non-zero also for an isotropic potential energy surface. Second, this 
is confirmed by the good quality of pressure broadening coefficients obtained by approximate 
semi-classical methods, where the impinging trajectory of the perturber is totally decoupled 
from the tensorial nature of the molecule/molecule electromagnetic interaction [3TH33] . In 
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retrospect, we understand that approximating the full Eq.(J3J) by the purely inelastic terms 
of Eq.Q may be correct in certain cases, but this should at least be carefully checked. An 
analogous point has been made earlier, see [23]. 

We show in figure [3] the present results as compared to estimates based on the inelastic 
sum, Eq. (4), but truncated to the first two terms, as used by DDP10 (and corrected for an 
error of a factor of 2 for the Dubernet et al. data). We see that there is a strong disagreement 
between the two computational schemes. In particular, data computed from the inelastic 
rates of Dubernet et al. [IS], which are based on the same PES as the present work, are 
significantly lower than the present rigorous calculations. The observed differences are to be 
attributed to the approximation in applying the Random Phase Approximation (neglect of 
the elastic scattering interference term) to Eq. (|4]). Furthermore, DDP10 made the further 
simplification of replacing the averaging procedure of Eq. |2jby using rate coefficients divided 
by the average thermal velocity. While this should not change the trend of the cr p b (T), it 
adds some further imprecision. For T > 20 K, the DDP10 points were extracted, via the 
same procedure, from the older values of Phillips et al. |34j based on a less accurate PES. 
We note that the same procedure was also applied by DDP10 to the case of H2O— He where 
it was found to be quite accurate, suggesting a very different structure of the scattering 
amplitudes. This is actually not surprising since i) the H^O-He and H2O-H2 PES are very 
different (see, e.g., [3j,ll]) and ii) additional coupling terms are introduced by the rotational 
angular momentum of H 2 . 

The low temperature range of the theory vs. experiment comparison, T < 80 K, is more 
problematic. Experiments show a dramatic decrease of cr p ' b (T) as T goes below about 80 K, 
for both transitions examined here and also for higher transitions. An analogous, even if 
less pronounced effect was also found for a p (T) in [T21 [13], for HCN-He and 13 CO-He 
respectively, albeit at much lower T, around 5K. No definite explanation may be found in 
these preliminary calculations, but it must be noted that disagreement between computa- 
tions and experiments arise at energies where the pressure broadening and inelastic cross 
sections enter in a regime where narrow resonances become prominent, see e.g. [28J. If res- 
onances are sufficiently narrow, that is, if the complex H2O — H2 is sufficiently long-lived, 
the impact approximation may be no longer valid. According to the density in the DDP010 
experiment (n ~ 10 20 cm -3 , B. Drouin, private communication), and with a cross-section 
of about a ph '(E = 80 K) ~ 50A 2 , an average speed of v = J2E/fj, ~ 1 km/sec, and a 
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FIG. 3: Pressure broadening cross-sections for the transition at 556 GHz, ortho-H 2 0, lio lot- 
Blue, red and black continuous lines, present theories, as in figure [TJ Open symbols, summation of 
inelastic scattering cross sections as in DDP10, circles from Dubernet et al., 2002, squares, from 
Green et al., 1996, with color scheme identical to the solid lines. Green triangles, experimental 
values taken from DDP10. 

resonance width of about T ~ lcm -1 [15] . we have that the interaction volume U becomes 
comparable to the inverse density [19] : 

U = (hva p - h -/T) ~ 1.6 1(T 20 cm 3 ~ l/n ~ 1(T 20 cm 3 , (5) 

invalidating the impact approximation. It is also possible that, for yet unknown reasons, 
there is a dramatic conversion from ortho-H 2 to para-H 2 at low temperatures, through 
some paramagnetic impurities in the walls of the cell [33], even if there is no experimental 
hint that indicates this explanation at the moment. Further insight is obviously needed in 
those regimes. Measurements with para-H 2 (j = 0) would be particularly valuable both for 
comparison with theory and for applications to cold interstellar clouds where H 2 is mostly 
in its para form [36] . 

In this Communication, we have shown using full quantum scattering methods combined 
with state-of-art PES, that very good agreement is found between experimental pressure 
broadening and theoretical calculations at temperatures where the impact approximation 
is clearly valid, i.e. T > 80 K. We have found that the rigorous theory of Baranger is 
needed to make valuable predictions for the present system and that simple approximations 
based on the manipulation of inelastic rates or cross sections are unreliable. We stress in 
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particular that the sharp drop observed in the pressure broadening measurements below 
~ 80 K, and not predicted by the present calculations, does not cast doubt on the inelastic 
rates computed by Dubernet et al. [5J [TBI ES] since a rigorous quantum theory of broadening 
beyond the impact approximation seems necessary in this regime, as discussed above. As a 
result, to our opinion, the water-hydrogen V08 PES is once more successfully tested, with 
a special emphasis on the mid- to long-range region of the PES and the isotropic part. We 
thus complement here various tests performed recently on the water-hydrogen system, such 
as differential measurements jl], and molecular beam scattering experiments [5], which have 
all so far confirmed the high accuracy of the V08 PES. Another series of experiments, now 
underway, would aim at the spectroscopy of the bound H2O-H2 van der Waals molecule. 
With all these tests completed in the near future, the V08 PES will be thoroughly tested 
and extremely good confidence should be gained for all astrophysical applications. 
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